Hadron Formation and the Phase Diagram of QCD Matter 
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We investigate hadronic species freeze-out systematics in A+A collisions at low SPS energies, 
corresponding to a baryochemical potential above 300 MeV, analyzing NA49 hadron production 
data in the framework of the statistical hadronization model, and in the UrQMD hadronic transport 
model. Observing no deviation from universal grand canonical hadro-chemical equilibrium freeze- 
out, we argue that the observed hadronic freeze-out points should universally signal the boundary 
line of the hadronic phase in the QCD matter phase diagram. 
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The phase diagram of strongly interacting mat- 
ter represents the most challenging open problem of 
quantum chromodynamics (QCD). Its key- feature, the 
confinement-deconfinement line between hadrons and 
partons interpolates between the regions in the [T, fis] 
plane where the temperature T is high and the bary- 
ochemical potential (1b (the net baryon density) low, to 
the domain of low T and high /is- Lattice QCD stud- 
ies at finite temperature and /j,b = have described the 
phase transformation as a rapid crossover in the vicin- 
ity of T c = 170 MeV [l| which coincides with Hage- 
dorns former estimate of the limiting temperature Th, 
expected for matter consisting of hadrons and their reso- 
nances Q. Toward higher /is, the position of the parton- 
hadron coexistence line has been estimated recently by 
extrapolations of lattice QCD to finite net baryon den- 
sity [H, 0, S], and by former considerations of the QCD 
process of chiral symmetry restoration, occurring at the 
border of hadronic matter 0, @, Q which is both confined, 
and exhibits strongly broken chiral symmetry. These two 
properties disappear at the QCD phase boundary. 

This phase boundary is traversed in relativistic 
nucleus-nucleus (A + A) collisions at sufficiently high yfs 
(above about 5-10 GeV). It is the goal of such studies to 
elucidate the properties of QCD matter as sampled along 
the dynamical trajectory of the collisional "fireball" . The 
dynamics in such collisions can be sketched 0, [ljj, [H| by 
an initial stage of matter heating and compression by pro- 
cesses occurring at the microscopic level (such as breakup 
of the parton structure functions), to be followed by a 
slowdown phase broadly analogous to a classical turn- 
ing point, during which partonic equilibrium should be 
approached and the system first appears on the [T, fig] 
plane of QCD matter, in various domains depending on 
A and \fs. The dynamical trajectory enters the phase 
diagram well above T c . A subsequent expansion phase 
transports the contained matter across the confinement 
and chiral symmetry breaking phase border(s), giving 
rise to hadronization, and hadronic "freeze-out" to a high 
T hadron-resonance system. A last stage, of essentially 
elastic and resonant rescattering, then leads to a final de- 
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FIG. 1: The phase diagram of QCD matter in the grand 
canonical variables of temperature T vs. baryochemicalpo- 
tential \ib ■ Also shown is a conjectured critical point [3, 3 IH] 
and the hadronic freeze-out points [l2l [l4l. HH [To] (see text). 



coupling from strong interaction. This so-called "kinetic 
freeze-out" occurs at lower temperatures, of about 100 
MeV, far below the phase boundary line. 

In this letter we shall focus on QCD confinement via 
hadronization, and its consequences observable in A + A 
collisions. The hadronization process results in an ob- 
servable ensemble of created hadrons and hadronic reso- 
nances, with a yield distribution over the various species 
that " freezes-out" in the immediate vicinity of the phase 
transformation [12, [H, [3, [HI, [IB]. Most remarkably, 
the systematic study of these measured yield distribu- 
tions, in terms of the statistical hadronization model, has 
revealed a universal "hadro-chemical" equilibrium, the 
yield distributions resembling grand canonical Gibbs en- 
sembles of hadrons and resonances [171, [l8|, [19, 20] . The 
derived parameters of the population, T and /j,b, vary 
monotonously with y/s. This model also succeeds in de- 
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scribing hadron production multiplicities in elementary 
e+e~ collisions !19j, suggesting that statistical equilib- 
rium production is a gen eral property of the hadroniza- 
tion process itself [H, [H H3, El, [H ill, HI • The result- 
ing idea - that the hadro-chemical freeze-out data in A+A 
collisions at various thus should result in a succession 
of points that universally locate the QCD hadronization 
transition in the [T, //g] phase diagram - will be pursued 
below. 

Fig. Q] illustrates our present view [3, [2j| H|| of the 
QCD phase diagram. A critical point on the deconfinc- 
ment line marks the end of the cross-over domain [l| 
at lower /ir, the transition becoming first order toward 
higher \xb @» 0> 0] ■ There, we also sketch a further phase 
transition line which defines the border between hadron- 
resonance matter and the quarkyonium phase recently 
postulated by McLerran and Pisarski [27|]. It is sug- 
gested to represent the large number of color (N c ) limit 
of QCD below T = T c , consisting of massive quarks plus 
thermal excitations of glueballs and mesons. Whereas 
the deconfinement line at low /is, and an estimate of 
the critical point position are adopted from recent lattice 
calculations 0], the border line of the quarkyonic phase 
is a mere guess. Following a suggestion by McLerran 
[111 we have sketched it such that it intersects the points 
corresponding to hadron-resonance freeze-out at various 
y/s 7 as determined from a statistical model analysis to be 
described below. We thus imply that hadron-resonance 
freeze-out occurs, at high /is, at the phase boundary of 
a new state of QCD matter, interpolating between the 
lattice confinement line and hadronization - such as the 
conjectured quarkyonic phase. 

This choice is in line with our above ascertation that 
hadronic freeze-out universally coincides with hadroniza- 
tion. This was first argued pj, Q [H [H, H, H only 
for low \ib where we see, in Fig. [TJ that the freeze-out 
points merge closely with the lattice phase boundary of 
the quark-gluon plasma (QGP). However, the freeze-out 
points fall well below this boundary at /is > 300 MeV, 
leading to the question of an intervening state of matter 
that could cool down the system toward hadronic freeze- 
out. 

It is the aim of this letter to show, first, that the grand 
canonical statistical model gives an excellent description 
of the hadron multiplicities obtained by NA49 [1^, [30] , in 
central Pb + Pb collisions at 40 and 30 AGeV (yfs = 8.7 
and 7.6GeV) where [Ib > 300 MeV. There is no deviation 
from a synchronous freeze-out in grand canonical species 
equilibrium. This indicates a freeze-out driven by a phase 
transition to the hadron-resonance state - as conjectured 
above. Second, we shall show that, in particular, the 
system state above the freeze-out points and below the 
confinement line can not merely be a more dense, hotter 
hadron-resonance matter as described, e.g., in transport 
models such as UrQMD [3l| because such a state would 
decouple to frozen-out hadron yields in a sequential pat- 
tern, in order of species specific inelastic mean free path. 
We thus seem to require a further non-hadronic phase 
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FIG. 2: Hadron production average multiplicities observed 
[HHj in central Pb+Pb collisions at AO A and 30A GeV (s/s 
= 8.7 and 7.6 GeV), vs. predictions of the grand canonical 
statistical model. 



above the hadronic freeze-out points, which decays into 
hadrons in a phase transition. The quarkyonic phase 
would provide for this. 

Turning to analysis of total A-k hadron production mul- 
tiplicities reported by the SPS experiment NA49 [lij [3(| , 
for central Pb+Pb collisions at 40 and 30^4Gey, we show 
in Fig. [2] the results of a statistical grand canonical en- 
semble model calculation. The formal and practical de- 
tails of this analysis are described in ref. [2(|. We note, 
in particular, that the so-called "strangeness undersatu- 
ration factor" 7 S < 1 is employed here, as a necessary 
ingredient of the fit procedure of total 47r multiplicities. 
It has been recently demonstrated for A+A collisions at 
SPS and RHIC [32|, [33| that this fudge factor may stem 
from surface effects: a significant fraction of the partici- 
pating nucleons from target and projectile nuclei interact 



3 



only once due to the dilute surface of nuclear density pro- 
files, thus experiencing canonical strangeness suppression 
characteristic of elementary p + p collisions (core-corona 
model) [33| . They do not participate in the bulk, grand 
canonical fireball. However, we can not make use of 
this quantitative model here as the elementary nucleon- 
nucleon corona cross sections are only poorly known at 
such low energies. 

The statistical model fits for the two energies illus- 
trated in Fig. O are quite satisfactory, with resulting 
[T, mb] = [148, 387] and [145, 434] and 7s = [0.86, 0.79] 
respectively. Their quality is equal to that encountered 
at higher e nerg ies, toward top SPS and RHIC energies 
[l7l \l§L [l9l |20| ]. This holds, in particular, for the pre- 
diction of hyperon yields which, thus, appear to emerge 
from freeze-out at the respective common bulk hadron 
temperatures. These temperatures fall well below the 
lattice QCD estimates of T c (see Fig. [1]), at these high 
values of but the system arrives at freeze-out in a 
global equilibrium of species, with no indications of a 
sequential freeze-out. 

The two resulting freeze-out points are included in 
Fig. [TJ The remaining points are taken from a former, 
similar analysis [34] . Our results for the chemical freeze- 
out temperature and baryon-chemical potential are in 
fairly good agreement with those presented in [l& l3a. |36|] . 

We conclude that the equilibrium distribution of 
hadron-resonance species at high \ib freeze-out reflects 
a system expansion trajectory from T c to Tf reeze — ou t, 
through a state of fireball matter that can cool down the 
species distribution while maintaining equilibrium. We 
shall show, next, that this state can not be, merely, a 
hadron-resonance population undergoing binary inelas- 
tic rescattering at the microscopic level, at density and 
temperature approaching T c . Such a system can not cool 
down the species distribution to a uniform equilibrium 
population, freezing-out at T well below Tc- 

This property of a dense hadron-resonance gas can 
be studied with the microscopic hadron-resonance trans- 
port model UrQMD [3l[. The absence of cooling the 
species distributions established at an initial phase tran- 
sition at T c (where the hadron-resonance population is 
generated via the Cooper- Frye formalism [371]). is illus- 
trated in Fig. [31 In the top panel we employ results 
of a former UrQMD calculation by Dumitru and Bass 
[38| , for central Au + Au collisions at top RHIC energy 
{yfs = 200 GeV). The hadronic yield distribution de- 
rived from grand canonical Cooper-Frye hadronization 
at an assumed Tc — 170 MeV is, on the one hand, 
shown as if it were emitted into vacuum, at this stage, 
without further interaction. Alternatively, a UrQMD bi- 
nary hadron collision expansion phase is attached, as an 
" afterburner" . Within exceptions concerning some in- 
evitable final state rearrangement of baryon and anti- 
baryon yields, the overall primordial yield distribution is 
preserved, frozen-out throughout the hadronic expansion 
phase. 

This analysis refers to y/s = 200 GeV where \ib is 
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FIG. 3: Top panel: Hadronic species multiplicity distribu- 
tions in central Au+Au collisions at y/s = 200 GeV result- 
ing, alternatively, from a model of direct freeze-out into vac- 
uum from hadronization at T = 170 MeV, and from an at- 
tached UrQMD hadronic expansion model as "afterburner" 
[3Sl ]. Bottom panel: Similar study of multiplicities in Pb+Pb 
central collisions at yfs = 8.7 GeV, with hadronization at 160 
MeV M. 



small (of order 20 MeV 18}). The bottompart of Fig. [3] 
illustrates results of a similar procedure [39| , carried out 
at y/s = 8.7 GeV, fi B = 380 MeV: the region of inter- 
est to the present study. The Cooper-Frye hadronization 
mechanism is employed, according to an average criti- 
cal energy density criterion, with resulting hadronization 
temperatures of 155-160 MeV. Again, the bulk hadronic 
species yields (e.g. ir, K, A and even E and Q) survive 
the hadronic expansion phase, as modeled by UrQMD, 
essentially unchanged. However, the anti-hyperon yields 
exhibit significant annihilation, in this case. This re- 
sults in a distortion of the global equilibrium distribution, 
which seems to be absent in the statistical model analy- 
sis illustrated in Fig. [2j We conclude, for now, that the 
mode of expansion modeled by UrQMD can apparently 
not cool down the system to a new equilibrium state, at a 
lower temperature, but seems to merely lead away from 
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an initial equilibrium. We leave these observations to 
the forthcoming, systematic investigation of the UrQMD 
model predictions. 

At low \ib the statistical model freeze-out points re- 
cover the parton-hadron phase boundary at T c (/j-b), as 
predicted by lattice QCD (see Fig. []}. The primor- 
dial hadron-resonance populations emerging at the QGP 
phase boundary thus are, apparently, not cooled down to 
below T c . This observation supports the hypothesis of a 
phase transformation from a QGP to a hadron-resonance 
phase [40], occurring at fis < 200 MeV. Conversely, 
however, the statistical model analysis at high (is de- 
termines a hadronization equilibrium temperature well 
below the lattice QCD estimate of T c , as we have seen 
in Figs. [T]and [2] Our above considerations show that an 
intervening hadron-resonance state of system expansion 
can not explain these observations. It can not cool down 
the bulk species distributions to below T c . A different 
state of QCD matter thus seems to be required, to shift 
the process of hadronization toward T(freeze-out) < T c 
at high /is, via a further phase transition. 

In summary, we have revisited the hadronic species 
freeze-out data at low SPS energies where \jlb ~ 
400 MeV, comparing them to the statistical model and 
to the UrQMD hadron-resonance transport model. The 
hadronic yield distributions do not exhibit aspects of a 
sequential freeze-out. On the contrary, freeze-out occurs 
in a global equilibrium grand canonical ensemble, similar 



to the situation encountered at [is — ► 0. We have argued 
that this can not be the consequence of a mere expansive 
dilution as it occurs in UrQMD. 

It has been pointed out 0, HI G3, [H, [H, 53 that 
hadronic freeze-out should coincide with hadronization 
of a quark-gluon phase, at low \ib- The universality of 
statistical equilibrium among the hadronic species sug- 
gests that the hadrons produced in A + A collisions uni- 
versally freeze out as the consequence of a hadronization 
phase transition, freeze-out thus always coinciding with 
hadronization. This implies that the observed freeze-out 
temperatures at high /i# also coincide with the appro- 
priate Hagedorn limiting hadronic temperature [42 1. At 
high fiB the phase boundary T c predicted by lattice QCD 
falls well above the freeze-out points as shown in Fig. [T] 
One thus expects a further, intervening phase of QCD 
which would postpone the hadronization transition, such 
as the quarkyonium phase predicted by Mclerran and 
Pisarski [27|,[28j]. 
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